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Following an earlier proposal by Choudhuri (2003) for the 
origin of twist in the magnetic fields of solar active regions, we 
model the penetration of a wrapped up background poloidal 
field into a toroidal magnetic flux tube rising through the so¬ 
lar convective zone. The rise of the straight, cylindrical flux 
tube is followed by numerically solving the induction equation 
in a comoving Lagrangian frame, while an external poloidal 
magnetic field is assumed to be radially advected onto the tube 
with a speed corresponding to the rise velocity. One predic¬ 
tion of our model is the existence of a ring of reverse current 
helicity on the periphery of active regions. On the other hand, 
the amplitude of the resulting twist depends sensitively on the 
assumed structure (diffuse vs. concentrated/intermittent) of the 
active region magnetic field right before its emergence, and on 
the assumed vertical profile of the poloidal field. Nevertheless, 
in the model with the most plausible choice of assumptions a 
mean twist comparable to the observations results. Our results 
indicate that the contribution of this mechanism to the twist can 
be quite significant, and under favourable circumstances it can 
potentially account for most of the current helicity observed in 
active regions. 
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Abstract. 


1. Introduction 


The magnetic field of a typical sunspot is approximately 
vertical below the photosphere and then spreads radially at 
the photospheric level. However, even visual observations 
show that the penumbral stru ctures of su nspot s are often 
twisted , as first pointed out by Hale ill927ll and iRichards^ 
lll94ll) . Vector magnetogram measurements in the past decade 
have established that sunspot magnetic fields have helical 
structures, with a higher occurrence of negative helicitv in 
the n o rthern hemisphere irSeehaferlll99f)t IPevtsov et al.lll99^ 


I20011 Bao & Zhang Il99^ . One possible theoret ical expla¬ 
nation of the observed helicity was proposed by IChoudhuril 
( I 2 OO 3 I) . who suggested that the poloidal flux in the solar 
conve ction zone (SCZ) gets wrapped around a rising flux 
tube. IChoudhuri et al. I J2nn4l) carried out a numerical sim¬ 
ulation demonstrating that the observed hemispheric hand¬ 
edness rule can indeed be reprod uced by incorporating this 
mechanism in their dynamo model ( iNandy & Choudhuril2002t 

IChatteriee et alj2004tlChoudhuri et alj2005). 


If the magnetic flux in the rising flux tube is nearly frozen, 
then we expect that the poloidal flux collected by it during its 
rise through the SCZ would be confined in a narrow sheath 
at its outer periphery. In order to produce a twist in the flux 
tube, the poloidal field needs to diffuse from the sheath into 
the tube by turbulent diffusion. However, turbulent diffusion is 
strongly suppressed by the magnetic field in the tube. This non- 
lin ear diffusion process was studied in an untwisted flux tube 
by Petrovav & Moreno-Inserti^ ill997lf who concluded that a 
substantial amount of flux may be eroded away from a rising 
flux tube during the process of its rise through the SCZ. The 
mod el was subsequently successfully appl ied for sunspot de¬ 
cay llPetrovav & van Driel-Gesztelvil 1997ll . In the present pa¬ 
per we extend this model by including the poloidal component 
of the magnetic field (i.e. the field which gets wrapped around 
the flux tube) and study the evolution of the magnetic field in 


the rising flux tube, as it keeps collecting more poloidal flux 
during its rise and as turbulent diffusion keeps acting on it. 

The conclusions drawn from our model hinge on some as¬ 
sumptions, especially concerning the subsurface magnetic field 
structure in the last phases of the rise of the tube. Allowing 
the various parameters of the model to vary over reasonable 
ranges, we find that the poloidal field wrapped around the flux 
tube should be able to penetrate inside the flux tube if the mag¬ 
netic field falls below the equipartition value in the top layers 
of the convection zone. On the other hand, if some physical ef¬ 
fects keep the magnetic field well above the equipartition value, 
then the poloidal field may remain confined in a sheath instead 
of penetrating to the core of the flux tube. We expect that much 
more high-quality magnetogram data will be available in the 
future and more will be known as to how the current helic¬ 
ity in a sunspot and in an active region varies with radial dis¬ 
tance from the centre. This may enable detailed comparisons 
between theory and observations in future, with the possibil¬ 
ity of constraining various parameters in the theory. Our aim 
at present is only to set up the basic theoretical framework and 
study some exploratory solutions. 

Section 2 presents the mathematical formulation of our 
problem, while the numerical solutions are presented and dis¬ 
cussed in Section 3. Section 4 concludes the paper. 


2. Mathematical formulation 

We consider a straight, cylindrical, horizontal magnetic flux 
tube rising through the solar convective zone. As all variables 
in this model depend only on the radial distance from the tube 
axis and on time, we will study the wrapping of the large-scale 
poloidal field around the flux tube by considering a radially 
symmetric accretion of azimuthal field by the flux tube. A fur¬ 
ther complication is the expansion of the flux tube during its 
rise, due to the decrease of the external pressure. This expan¬ 
sion is assumed to be self-similar. Certainly self-similarity is 
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not a bad assumption. We expect the tube to expand in such a On substituting (6)-(l 1) into (4) and (5), a few steps of straight¬ 
way that the density inside the tube remains homogeneous. A forward algebra give 
self-similar expansion ensures that. 

2.1. Equations in the comoving Lagrangian frame 


dK 

dt 


d I dB' 


( 12 ) 


Suppose we formulate our problem in a frame of reference 
fixed with the centre of the rising flux tube. Substituting a radial 
expansion velocity 

V = VCr (1) 

and an axisymmetric twisted magnetic field 

B = B-Cj + (2) 


in the induction equation 

BB 

— = V X (v X B) - V X ( 77 V X B), (3) 

dt 

we get the following equations for the axial and poloidal fields 
in the flux tube; 


dB^ \ d , 


1 _B 

~r¥r 



( 4 ) 


BB^ B B 

^ ^ ^ Jr 


1 d 

Tl--^(.rB^) 
r or 


( 5 ) 


We first consider what happens in the interior of the flux 
tube as it expands and is subject to turbulent diffusion. We shall 
discuss later how accumulation of additional poloidal flux dur¬ 
ing its rise can be incorporated. The independent variables in 
the two equations (4) and (5) are r and t. Let us assume that the 
material inside the flux tube expands in a self-similar fashion 
and use the Lagrangian position coordinate ^ of a fluid element 
instead of r. The initial value of r at time t - 0 can be taken as 
the value of Then we can write 


f = F{ty, 


(6) 


where F{t) will have to be a monotonically decreasing func¬ 
tion of t for an expanding flux tube. The assumption of self¬ 
similarity implies that we have to use the same factor F{t) for 
all fluid elements inside the flux tube to go from the radial co¬ 
ordinate r to the Lagrangian coordinate In equations (4) and 
(5) we now want to transform from variables (r, t) to (^, t). By 
the chain rule of partial differentiation, it can easily be shown 
that 




The velocity is given by 



V - 


dt 




( 8 ) 

( 9 ) 


We also substitute 


B,(r,t)^B'^{^,t)FHt), ( 10 ) 

Byr,t)^B'^(^,t)F(t). ( 11 ) 




1 d 

'77 ^(^^P 


( 13 ) 


It is thus clear that B'(^, f) and B^(^, f) do not change with time 
if 77 = 0 , in accordance with what we expect under the condition 
of flux freezing. 

The equations (12) and (13) clearly should hold from the 
centre to the outer periphery of the flux tube if we assume it to 
expand self-similarly. 

Now we turn to the problem of how to incorporate in our 
equations the accretion of poloidal flux to the tube. As the flux 
tube rises, it collects more poloidal flux, which gets wrapped 
around it. If we are in the frame of the rising flux tube, it would 
seem that there is a flow of fluid from the upward direction 
bringing the poloidal flux. In the downward direction, the sur¬ 
rounding fluid flows away from the flux tube. However, the ten¬ 
sion of the poloidal flux makes sure that it gets wrapped around 
the flux tub e, after the reconn ection in the wake, as shown in 
Figure 4 of IChoudhuril ll20()3l) . In our one-dimensional model, 
we can approximately take account of this by assuming that the 
poloidal flux is brought uniformly from all directions by a ra¬ 
dial inward flow with velocity equal to the velocity with which 
the fluid is flowing from the upward direction. Let us, therefore, 
consider the nature of flow velocity from the upward direction. 
Finding the flow past a cylinder is a standard problem in incom¬ 
pressible fluid dynamics (which should hold under the subsonic 
conditions prevailing in our problem) and is discussed in many 
standard textbooks (see, for example, IChoudhurilll998i §4.7). 
To get the velocity in the upward direction, we basically have to 
substitute 0 = 0 in the expression (4.54) of IChoudhuril 1119981) 
giving the velocity potential and then differentiate it with re¬ 
spect to r. This gives 


v^-Uf, 1 


J 


(14) 


where Uft is the velocity of rise of the flux tube and ^ 1/2 is its 
radius, defined as the place where B^ falls to half its maximum 
value (at the centre of the flux tube). We assume that there is an 
isotropic radial inward flow given by equation (14) towards the 
flux tube from all directions in the region ^>^ 1 / 2 , whereas the 
flow is zero inside the flux tube (^ < ^ 1 / 2 ). Such a flow field 
certainly has a non-zero divergence and other bad properties. 
But it should capture the basic physics of poloidal flux advec- 
tion in a one-dimensional model. We now discuss how we And 
the advection of the poloidal flux by such a flow field. The term 
d(vBf)ldr in equation (5) gives the advection of the poloidal 
held by a radial flow. While deriving equation (13) from equa¬ 
tion (5), we got rid of this term by assuming the self-similar 
expansion which gives the velocity held (9). If there is an ad¬ 
ditional velocity held, then the advection term for that velocity 
would still persist. So we put an additional advection term in 
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(13) in our scaled variable, so that equation (13) becomes 


OB' 






(15) 


dt 

The factor F in the last term comes from equation (7). 


2.2. Input parameters: rise of the flux tube 

To understand the magnetic held evolution in the hux tube, we 
have to solve equations (12) and (15), starting from some initial 
conhguration and setting suitable boundary conditions at the 
outer periphery of our region of integration that would allow 
the free inward advection of the poloidal hux. We also need 
to specify F, rj and v. Of these, v has already been specihed 
by equation (14), although we need to know Uft to get v. We 
now discuss how we obtain F, Uft and rj in our model. We 
shall come to a discussion of the initial and boundary condi¬ 
tions later. 

As the hux tube rises through the SCZ, we denote its radial 
distance from the centre of the Sun by R. The hux tube begins 
from the bottom of SCZ at = Rh, where its radius is and 
the external density is When it rises to R where the external 
density is pe, its radius becomes rp. Since the density inside the 
Hux tube would be very nearly equal to the external density, 
mass conservation implies 


RhPe.off, = RPerjf 
From (6) it follows that 


p ^ I RPe 
rp \ RbPefi ' 


(16) 


Thus, to hnd F as a function of time, we need to hnd out how 
R changes as a function of time and we also require a model of 
the SCZ which will give us the value of pe at that value of R. 
We now discuss how we prescribe a model for SCZ and how 
we calculate the rise of the Hux tube through this SCZ, giving 
/? as a function of time. Since Uft appearing in equation (14) is 
essentially given by dRIdt, the velocity v as given by equation 
(14) also gets completely specihed once we know how R varies 
with t. It is also clear from (16) that F would have a specihc 
value at a certain depth R within the convection zone. The solid 
line in Figure 1 shows F as a function of R. 

How a horizontal magnetic Hux tube rises through SCZ can 
be studied in a fairly straightforwa rd fashion jMoreno-lnsertisI 
Il983t Choudhuri & Gilman Il987ll . We believe that active re¬ 
gions form by the buoyant rise of a part of a hux tube, while 
other parts remain anchored at the bottom of SCZ. Studies of 
the rise of such loops show that the upper p arts of the loops 
move very much like horizontal Hux tube s llMore no-lnserti^ 
' l986tlChoudh^[l989t5’Silva&Choudhurilll993ll . We hnd 


out how R varies with t by considering the rise of an ax- 
isymmetric Hux ring. The dynamics of such Hux ring s has 
been studied exhaustively bv ichoudhuri & Gilmanllfl98^ . The 
forces acting on such a hux ring are (i) magnetic buoyancy, 
(ii) the Coriolis force, (iii) magnetic tension, and (iv) the drag. 
The simulations match various aspects of observational data 


best if the initi al magnetic held at the bottom of SCZ is as 
strong as 10^ G (IChoudhuri & Gilma Jl987HChoudhurill989t 

IP’SHva & Choudhurill993l:lFan et al.ll993l) For such a strong 

initial magnetic held, the Coriolis force is un important and we 
neglect it. Also. lchoudhiiri & nilmanlill987l) made an estimate 
of magnetic tension compared to magnetic buoyancy. Although 
magnetic tension may be an appreciable fraction of the mag¬ 
netic buoyancy at the bottom of SCZ, it becomes negligible 
in the upper parts of SCZ. So we neglect magnetic tension 
also. The neglect of the Coriolis force and magnetic tension 
makes sure that the hux tube moves radially and the problem 
becomes one-dimensional. T he only terms in equation (4) of 
IChoudhuri & GilmanI lll987l) which should be of interest to us 
are 


, d^R , 

= -{nii - Me 



-jCoPerp 



(17) 


where gs is the gravity at the solar surface, whereas m,- = ^rr^^p, 
and lUe - J^r^pPe are respectively mass per unit length of the hux 
tube and the displaced huid. The dimensionless drag coefficient 
Co is f ound to have a value of around 0.4 in laboratory exper¬ 
iments (lGoldsteinlll938tlSchlichtinglll979l) . Dividing equation 
(17) by Inr^pPe, we get 



It may be noted that pi ^ p^ so that we treat them differently 
only when we have to consider the difference between them. 

The magnetic buoyancy factor (p^ - p,) /(2p(,) depends on 
the background model of SCZ that we use. If we assume 
the temperature gradient to be exactly equal to the adiabatic 
gradient, then we end up with a polytropic model for SCZ. 
Choudhuri & Gilman (1987) used a polytropic model with a 
choice of parameters which gave a close ht to more detailed 
models of S CZ. The model is de s cribed through Equations 
(10)-(12) of IChoudhuri & GilmanI (Il987h . with the values of 
parameters listed at the beginning of § III. We use this model 
of SCZ here. We also need to use some suitable thermal condi- 
tion to calculate the magneti c buoyancy factor (p^ - p,)/(2p(,). 
IChoudhuri & GilmanI lll987l) presented results for three ther¬ 
mal conditions. We present our calculations here for the sim¬ 
plest case of the Hux tube being in thermal equilibrium with 
the surrounding, although we have done some calculations for 
the other thermal conditions and found the results to be qualita¬ 
tively similar. For the case of thermal equilibrium, the magnetic 
buoyancy factor is given by 


Pe-Pi ^ Bl /^y2-r)/(r-i) 

2pe IhnpefiXTo) [Rbl 


(19) 


Knowing the temperature T at the position R from the SCZ 
model, we can use equation (19) to calculate the magnetic 
buoyancy factor of a hux tube at R, which started with initial 
magnetic held Bq. Calculating magnetic buoyancy in this way 
and using equation (16) to hnd rp at the position R, we can 
integrate equation (18) to hnd how R changes with time. From 
the variation of R with time, on making use of equation (16), 
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order. It may be noted that we use equation (22) only when 
rj, < H, which is the case for typical flux tubes rising through 
the bulk of SCZ except the uppermost layers. There we take 
t/o = t/oo- However, even the validity of equation (18) becomes 
questionable if rf, > H. 

Now equations (12) and (15) are a set of nonlinear flux- 
conserving equations that may be solved by an explicit two- 
step Lax-Wendroff scheme. The time step At obeys the stabil¬ 
ity condition At < mm[(A^)^Irjo, A^IUp], where U/t is the ve¬ 
locity of rise calculated from equation (18) and rjo is given by 
equation (22). We use a non-uniform but steady spatial grid 
with 1500 points which has a finer resolution of ~ 2.4 km 
for ^ < 3000 km and increasing successively by 2% thereafter 
upto ^ = 20,000 km. 

2.4. Initial and boundary conditions 


Fig. 1. Variation of the expansion factor F with distance R from 
the centre of the Sun during the rise of a horizontal flux tube in 
the convective zone (modelled as an adiabatic envelope). Solid; 
F calculated from eq. (16). Dotted and dashed-dotted; F sub¬ 
ject to the requirement B > Z?eq and B > 3Beq, respectively. 


we can And how F{t) varies with time (Fig. [Q. The velocity 
Uft in equation (14) is essentially the rise velocity of the flux 


tube, i.e. 



( 20 ) 


So, once we know how R changes with f, we can use equations 
(14) and (20) to obtain v that appears in equation (15). 


2.3. Input parameters: turbulent diffusivity 


As we already pointed out, we need to keep solving equa¬ 
tions (12) and (15), as the flux tube rises, to And out how 
the magnetic field evolves. Now that we know how to find 
F{t) and v at any time step, we only need to specify the 
turbulent diffusivity n. For this w e follow the approach of 
IPetrovav & Moreno-Inserti^ ll 19971) and take ij to be given by 
the expression 


70 

1 + |B/Beql" 


( 21 ) 


Now that we have described how all the various terms appear¬ 
ing in equations (12) and (15) are specified or can be obtained, 
we only have to discuss the initial and boundary conditions 
used to solve these equations. Suppose we want to consider 
the buoyant rise of a flux tube carrying initial flux 'T with an 
initial magnetic held Bq at the bottom of SCZ. The initial ra¬ 
dius of the flux tube is obviously ^/r = a/T'/ttBo. We take the 
initial condition that = Bq inside r < and B^ = 0 out¬ 
side, whereas Bf is initially taken to be zero everywhere. The 
integration region over which equations (12) and (15) are in¬ 
tegrated extends to rout which is typically taken at 10^/f. The 
solutions are not very sensitive as to what boundary conditions 
we use for B'. However, since the flux tube keeps on acquir¬ 
ing poloidal flux which must be advected inward through the 
boundary of the integration region, the boundary condition on 
B^ is quite important. We assume that the SCZ is filled with a 
uniform horizontal magnetic held of 1 G, which gets wrapped 
around the rising flux tube. In order to achieve this, we have to 
continuously advect B^ = 1 G through the outer boundary of 
the integration region. Looking at equation (11), we realize that 
the appropriate boundary condition at r = rout is going to be 

B; = l/B(f)G. (24) 

3. Results and discussion 


where B = y B^ -H B^ is the amplitude of the magnetic held and 
Beq is the equ ipartition rnagnetic held. We use the convection 
zone model of ll Inno et al.l a 198.51) to obtain Beq at different posi¬ 
tions R within SCZ. Most of our calculations are done by taking 
X = 2 (except in the cas e presented in Figure 3). We specify rjo 
exactly the same way as [Petrovav & Moreno-Insertis|(ll997l) . If 
FI is the pressure scale height, then we take 

70 = 700I (22) 

with 

7700 = 3 X lO^^cm^s"^ (23) 

IChatteriee et alJ ll2004t) found that solar dynamo models give 
the best fits with observations for a value of diffusion of this 


Flux t ube simulations llC houdhuri & Gilrnan|l987t lChoudhu^ 


Il989t ID’ Silva & ChoudhuriHl 993t iFan et al]ll9'^ have sug¬ 


gested that theory matches various aspects of observations best 
if the flux tubes start with initial fields of order 10^ G at the 
bottom of SCZ. We, therefore, carry out all our calculations by 
taking an initial magnetic held of 10^ G. The typical flux car¬ 
ried by a large sunspot is 10^^ Mx. So we present results for a 
flux tube with such flux, which implies that the initial radius at 
the bottom of SCZ is Rt, - 1.78 x 10^ km. We study the rise of 
the flux tube by numerically integrating equation (18) with the 
magnetic buoyancy given by equation (19). While the flux tube 
rises, we study the evolution of the magnetic held in the flux 
tube by solving equations (12) and (15) using a Lax-Wendroff 
scheme, with F{t) given by equation (16), v by equation (14) 
and 77 by equation ( 21 ). 
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Fig. 2. Plots of B', B'^ and Up as functions of ^ for a rising flux 
tube with x - 2 (case A). The different curves correspond to 
the profiles of these quantities at the following positions of the 
flux tube: O.IRq (thick solid), O.SSBq (solid), 0.9Bo (dashed), 
O. 95 B 0 (dotted), O. 98 B 0 (dash-dotted). The flux tube reaches 
these positions at times 0 days, 5.3 days, 6.6 days, 7.9 days 
and 8.7 days after the initial start. The values of (B/Bgq) at the 
centres of these flux tubes at these positions are 10, 1.72, 0.64, 
0.098, 0.015. 

Figures 121 to |6l show the magnetic fields of the rising flux 
tube at depths R = OJRq, O.SSRq, 0.9Rq, 0.95Rb, 0.9SRq. 
Profiles of B' as functions of ^ are shown in the top panels, 
whereas the middle panels show the profiles of B^. As these 
panels show the rescaled variables, the actual held strengths 
can be calculated from them using Figure [Q to read off F for 
each curve, and plugging that value into equations (10) and 
(11). Note that B^ tends to the value given by (24) for large 
Since this is much smaller than what B^ becomes near the 
flux tube boundary, it appears in Figures|2lto|3as if B^ is going 
to zero for large although this is not the case. 

The bottom panels provide plots of 

ap^(VxB)JB, (25) 

since this is the quantity that essentially all photospheric 
measurements of the current helicity actua lly determine 
( iLeka & Skumani(^ll999t iBurnette et al.ll2003) . It is easy to 
see that ap, which has the dimension of 1/length, is invariant to 
the rescaling of the radial coordinate in our flux tube. 

3.1. Case A: diffuse magnetic fieid near the surface 

It should be evident from Figure|2lthat diffusion does not play 
a significant role until the flux tube reaches O.85B0. The diffu¬ 
sion in the deeper layers turns out to be negligible due to two 
factors: (i) since (B/Bgq) is large, the quenching included in 
equation (21) is quite efficient; (ii) since (rf,IH) is small, rjo 
as given by equation (22) turns out to be small. As a result of 
the low diffusion, the toroidal field B' does not change much, 
whereas the poloidal flux B^ remains confined in a sheath at 


the outer periphery of the flux tube (since the low diffusion 
does not allow B^ to diffuse inward). After the flux tube has 
risen above O.85B0, the toroidal magnetic held inside the flux 
tube becomes comparable to the equipartition held Bgq. Then 
diffusion is able to affect the magnetic held much more, since rj 
given by equation (21) becomes much larger due to the dimin¬ 
ishing role of magnetic suppression of diffusion. As a result, 
B' starts diffusing while the flux tube rises through the top lay¬ 
ers of SCZ, whereas B^ penetrates to the core of the flux tube 
instead of remaining confined to a sheath at the other periphery. 

In the very top layers, however, the cross-section of the flux 
tube becomes enormous in our model, leading to small values 
of F(t) and again making diffusion less important, as can be 
seen from equations (12) and (15). This is why we And that the 
magnetic held profile has not evolved that much from Q.95Rq 
to O. 98 B 0 . It may be noted that the profile of B' no longer has 
a sharp edge after the flux tube rises above O. 9 B 0 , but spreads 
around. In the realistic situation, however, we expect the evolu¬ 
tion of B' in the immediate vicinity of the flux tube to be much 
more complicated. In the frame of the flux tube, it would appear 
that the surrounding fluid is flowing inward from the upward 
direction and this would prevent the spread of B' in the upward 
direction. On the other hand, the surrounding fluid moves away 
from the flux tube in the downward direction and would carry 
away B' with it, leading to a larger spread of B' in the down¬ 
ward direction. In a one-dimensional model, the best way of 
capturing the average behaviour of B' may be not to include any 
velocity held in the surrounding fluid, as we have done by not 
including any advection term in equation (12). However, an ad- 
vection term is included in the evolution equation (15) for B^. 
Because of the topology of magnetic held lines, the poloidal 
held cannot be advected freely in the down ward direction, as 
can be seen in Figure 4 of ichoudhuril ll2003h . After a reconnec¬ 
tion in the wake, a poloidal held line should remain wrapped 
around the flux tube due to its tension. Hence the behaviour of 
the poloidal held is best captured in a one-dimensional model 
by including a uniform inward flow from all the directions. 
To treat the evolution of the magnetic held more realistically, 
it would be necessary to go beyond one-dimensional models. 
Some of us are now involved in developing a two-dimensional 
model of this problem. However, the one-dimensional model 
should capture some of the basic effects correct to the right or¬ 
der of magnitude. 


3.2. Parameter dependence 


Since many parameters of the problem are not known well, 
one very important question is whether the results presented 
in Figure 12] are sufficiently generic or would change substan¬ 
tially on changing the various parameters. Figure |2| presents 
results obtained by assumi ng thermal equilibrium of the flux 
tube with the surroundings. Choudhuri & r.ilmanllll987h pre¬ 
sented results for two other thermal conditions: by assuming 
the interior of the flux tube as adiabatic with and without a su- 
peradiabatic gradient in the surrounding SCZ. We carried out 
some calculations with these two thermal conditions as well. 
There are virtually no changes till the flux tube rises beyond 
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Fig. 3. Same as Figure|3but with x - 5. The flux tube reaches 
positions O.O.IRq, 0.85Rq, 0.9Rq, 0.95Rq, 0.98Rq at times 0 
days, 5.2 days, 6.6 days, 7.9 days and 8.7 days after the initial 
start. The values of (B/Beq) at the centres of these flux tubes at 
these positions are 10, 1.74, 0.66, 0.095 and 0.018 respectively. 



Fig. 4. Same as Figure |2] except t/oo = 10'^ cm^/s. The flux 
tube reaches positions O.IRq, 0.85/?©, 0.9Rq, 0.95Rq, 0.9SRq 
at times 0 days, 5.2 days, 6.6 days, 7.9 days and 8.7 days after 
the initial start. The values of (Z?/Z?eq) at the centres of these 
flux tubes at these positions are 10, 1.72, 0.69, 0.14 and 0.03 
respectively. 


0.9Rq. After that, the flux tube rises faster in these two cases 
compared to the case of thermal equilibrium and diffusion has 
less time to act. Therefore, we find the effect of diffusion a lit¬ 
tle bit less when the flux tube reaches the topmost layers of 
SCZ. Figure 13 presents results obtained by using a quenching 
index x = 5 instead of = 2 used in Figure |3 whereas all the 
other things remain the same as in Figure |2] Again the results 
are not qualitatively different. We also study what happens if 
the diffusivity is made smaller compared to what is prescribed 
in (23). Figure|3presents results obtained with 7700 = 10^^ cm^ 
s ', whereas all the other things remain the same as in Figure|3 
As expected, we find the diffusion somewhat less. Otherwise, 
results are not qualitatively different. We thus conclude that the 
results presented in Figure |2l are sufficiently generic for a rea¬ 
sonable range of parameters, as long as we hold on to our basic 
model of the horizontal flux tube rise. 

3.3. Case B: concentrated magnetic field near the 
surface 

One of the unsatisfactory aspects of the horizontal flux tube rise 
model is that the magnetic field falls to very small values in 
the top layers of SCZ when the flux tube expands enormously 
by moving to a low-density region. In the case presented in 
Figure |3 our model predicts that the magnetic field inside the 
flux tube when it reaches 0.98 Rq is about 400 G—an order of 
magnitude smaller than a typical sunspot field. Since such a 
field is much weaker than the equipartition field, it should be 
clear from equation (21) that the magnetic quenching of dif¬ 
fusion becomes completely negligible and diffusion is able to 
act without being inhibited by the presence of the magnetic 
field. The presence of 3000 G magnetic fields in sunspots is a 
compelling proof that magnetic fields may never fall to such 
low values; in fact, at least in photospheric layers, they remain 


well above the equipartition value (the magnetic field inside 
a typical sunspot being about thrice the eq uipartition field) . 
The non-axisymmetric flux tube simulation s llChoudhurill989t 
iD’Silva & Choudhurill993UFan et al.ll993h show that the fluid 



0.5 

0 

'e-0.5 

1 ^ 1.5 

s ‘^-2 



Fig. 5. Plots of B', B'^ and Op as functions of ^ for a rising 
flux tube. The field inside the tube is not allowed to decrease 
below Beq (case Bl). The different curves correspond to the 
profiles of these quantities at the following positions of the 
flux tube; O.IRq (thick solid), 0.85Bo (solid), 0.9Rq (dashed), 
O. 95 B 0 (dotted), O. 98 B 0 (dash-dotted). The flux tube reaches 
these positions at times 0 days, 5.2 days, 6.6 days, 7.7 days and 
8.2 days after the initial start. The values of (B/Beq) at the cen¬ 
tres of these flux tubes at these positions are 10, 1.72, 1.0, 1.0, 
1.0 respectively. 
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to a value sB^q at some depth (i being a numerical factor of the 
order of unity). We assume that the magnetic held inside the 
hux tube remains sB^q in the higher layers as it rises further. 
Here Bgq is the local equipartition value at the particular depth. 
If this is the case, then magnetic buoyancy would be given by 

^ ( 26 ) 

2pe I6npe 

instead of equation (19). While we calculate the rise of the Hux 
tube by using this expression of magnetic buoyancy, we cannot 
allow the cross-section to expand indehnitely if the magnetic 
Held has to remain sB^q. Instead of equation (16), we calculate 
F{t) by using the relation 

( 27 ) 


Fig. 6. Plots of B[, B'^ and ap as functions of ^ for a rising flux 
tube. The Held inside the flux tube is not allowed to decrease 
below 3Beq at any height (case B3). The different curves corre¬ 
spond to the profiles of these quantities at the following posi¬ 
tions of the flux tube: 0.1 Rq (thick solid), 0.85Rq (solid), 0.9Rq 
(dashed), 0.95Rq (dotted), 0.98Rq (dash-dotted). The flux tube 
reaches these positions at times 0 days, 5.2 days, 6.2 days, 6.8 
days and 7.0 days after the initial start. The values of (BfB^q) 
at the centres of these flux tubes at these positions are 10, 3, 3, 
3, 3 respectively. 


drains from the tops of rising magnetic loops making the mag¬ 
netic fiel d there strong er. Additionally, effects like convective 
collapse (l■SteineT^l2003l) can be operative near the solar s urface 
to enhance the magnetic field. iLongcooe & Choiidhiiril il2002t) 
have argued that flux tubes get distorted by convective turbu¬ 
lence in the top layers of SCZ, leading to the observed scatter 
of tilt angles around what would be expected from Joy’s law. 
Such buffeting of flux tubes by turbulence can also cause an 
enhancement of magnetic field by stretching. 

Most of the flux rise simulations are based on the thin 
flux tube approximation, which should be valid during the 
rise of flux tubes to about Q.9Rq during which the magnetic 
field remains sufficiently strong to stay relatively unaffected by 
the surrounding turbulence. So we can presumably trust the 
flux rise simulations through the deeper layers of convection. 
However, our understanding of what happens during the rise 
of flux tubes through the top layers of SCZ is extremely poor. 
Techniques of local helioseismology, such as time-distance 
seismology, have now evolved to a stage when the direct study 
of subsur face structures in emerging active regions ha s become 
possible jKosovichev et al.koOfltlHiighes et al.l200.5h . There is 
thus hope that this issue may be resolved in the not too distant 
future. In any case, the observation of sunspot magnetic fields 
indicates that there may be effects which prevent the magnetic 
field from falling to values lower than the equipartition field 
even in the top layers of SCZ. 

We now present some calculations by artificially not al¬ 
lowing the magnetic field to fall below the equipartition value. 
Suppose the magnetic field in the interior of the flux tube falls 


(dashed and dash-dotted lines in Figure^. By calculating F(t) 
in this way, we solve equations (12) and (15) to find the evolu¬ 
tion of the magnetic field. It may be noted that in the expression 
of drag in equation (18) also, we have to use rp = ^ftlF{t) with 
F(t) given by equation (27). 

Results for i = 1 (case Bl) and i = 3 (case B3) are shown 
in Figures|3and|^respectively. As in Figures|2j0 we plot B', 
B'^ and ap as functions of ^ at depths Q.IRq, 0.85/?©, Q.9Rq, 
0.95^0 and 0.98/?©. The times taken to reach these depths are 
given in the figure captions. Comparing with the times given 
in the caption of Figure |2] it will be seen that the flux tubes 
have risen faster through the top layers of SCZ, since magnetic 
buoyancy has remained stronger. The diffusion remains signifi¬ 
cantly quenched if the magnetic field stays higher than Bgq and 
has also less time to act because the flux tube rises faster. As a 
result, we see that the effect of diffusion is somewhat less for 
case Bl and drastically less for case B3 compared to case A. 
We see in Figure |6l that B' has not diffused much and B^ has 
remained confined in a narrow sheath at the boundary of the 
flux tube, being unable to diffuse inward. 


3 . 4 . Current helicity 


Observations of the current helicity parameter Op indicate that 
its typical average value in an active region is on the order of 
10“^ m~* l| PeY^?Py et al]ll995t Ivan Driel-Gesztelvi et al.ll2003l: 


iBurnette et al.ll2004 ). Inspecting the lower panels of Figures |2 

to 13 one finds that the typical value of Op in the internal parts 
of the flux tube is of order ~ 10“*m ' at a depth of 0.85/?© 
in all the cases studied. However, as the flux reaches the so¬ 
lar surface, in all the cases except the case B3 presented in 
Figure 6, the B^ component spreads out due to diffusion and 
its gradient becomes smaller, reducing ap by about one order 
of magnitude. Only if the magnetic field inside the flux tube 
remains stronger than the equipartition field (the case B3 rep¬ 
resented in Figure 6), the B^ component is unable to diffuse 
inside so that its gradient remains strong and ap is of order 
~ 10“* m“' even near the surface. This suggests that our case 
B3 may be closest to reality, i.e. during the rise of the flux tube 
from 0.9B© to 0.98B© effective flux concentration processes 
are at work, keeping the field strength at a value somewhat 
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above the equipartition level. This is consistent with the notion 
that the main flux concentration effect at work here is turbu- 
lent concentration (i.e. flux expulsion by the turbulent eddies: 
IProctor & Weissj [l982|l. while convective collapse and thermal 
relaxation ( Steineil2003ll are restricted to the shallowest layers 
above 0.98 Rq. 

Note, however, that in the calculations presented here the 
amplitude and sign of the poloidal field was assumed not to de¬ 
pend on depth. For alternative assumptions, significantly differ¬ 
ent current helicities may result, so the above conclusion should 
be treated with proper reservation. Details of the radial depen¬ 
dence of the poloidal field strength may strongly depend on the 
dynamo model. 

A more robust feature of the current helicity distributions, 
present in all the lower panels of our plots, is the presence of 
a ring around the tube with a current helicity of the opposite 
sense. This is clearly the consequence of the fact that on the 
outer side of the accreted sheath the radial gradient of the az¬ 
imuthal field, and thus the axial current, is negative. This is 
an inevitable corollary of the present mechanism of producing 
twist in active regions. A rather strong prediction of this model 
is, therefore, that a ring of reverse current helicity should be 
observed on the periphery of active regions, somewhere near 
the edge of the plage. 



Fig. 7. Flux loss as a function of height R inside the convection 
zone for the four cases presented in Figure|2l(solid),|3dashed), 
13 dotted), l^dash-dotted). 


calculations are repeated for the flux tube rising adiabatically 
through the convection zone, then the flux losses are somewhat 
less. For example, in the case B3, the flux loss is then 41.2% 
instead of 46.5% during the isothermal rise. 


3.5. Flux loss from the rising tube 

Petrovav & Moreno-Inserti^ ill9971) came to the conclusion that 
a considerable amount of magnetic flux is lost during the rise of 
a flux tube through the SCZ (see their Fig. 10). This conclusion 
was based on an approximate expression of inward velocity of 
turbulent erosion and involved various uncertainties. One im¬ 
portant question is whether our more careful calculations show 
similar flux losses. To make an estimate of the flux loss, we 
keep calculating the flux 

^0 = 27t (28) 

within the initial Lagrangian radius of the flux tube. Figure^! 
gives plots as i/tq as a function of the position in the SCZ for 
the four cases presented in Figures El El El and 13 We find 
that the flux loss is much less if the magnetic field is not al¬ 
lowed to fall below Bgq- Even in the case when the magnetic 
field is allowed to fall to rather low values (the solid line in 
FigurefTt, the flux loss is nowhere a s substantial as inferred by 
Petrovav & Moreno-lnserti^ lll997l) . The main reason for this 
discrepancy is that we do not make use of the assumption of 
a continuous “re-initialization” of the decay, as it was done in 
that paper. The basis of that assumption was that the flow of 
external fluid relative to the tube would instantly remove all 
magnetic flux lost from the tube. However, in the presence of 
an azimuthal field component, field line topology is expected 
to inhibit such flux removal. We thus conclude that flux loss 
during the rise of a flux tube is muc h less significant than what 
IPetrovav & Moreno-lnserti^ll 19971) found in their calculations. 
This provides a justification of the flux rise calculations based 
on the thin flux tube equation, where it is assumed that the mag¬ 
netic field is frozen during the rise of the flux tube. If the above 


4. Conclusion 

Alternative mechanisms for the origin of twist in active re¬ 
gion magnetic fields include buffeting of the rising flux tubes 
by helical turbulent motions jTongcone et al.lll998|i: the effect 
of Co riolis force on flows in rising flux loops ilFan & Gona 
l2flflf)l) : differential rotation (|De\^igl2QQ-(]l)j and h elicity gener¬ 
ation by the solar dynamo jSeehafer et alJE003l) . The relative 
importance of these processes is currently a subject of debate 
dHolder et al.ll2004l) . It is likely that the accretion of poloidal 
fields during the rise of a flux tube is just one contribution to the 
development of twist. Its importance may also be reduced by 
3D effects: considering the rise of a finite flux loop instead of 
an infinite horizontal tube, the possibility exists for the poloidal 
field to “open up”, giving way to the rising loop with less flux 
being wrapped around it. It is left for later multidimensional 
analyses of this problem to determine the importance of any 
such reduction. In any case, the results presented above indicate 
that the contribution of poloidal field accretion to the develop¬ 
ment of twist can be quite significant, and under favourable 
circumstances it can potentially account for most of the current 
helicity observed in active regions. 

In our calculations we found that while the flux tube rises to 
a depth of about 0.85Ro, the effect of diffusion is small and the 
poloidal field remains confined in a narrow sheath at the periph¬ 
ery of the flux tube. Afterwards, if the magnetic field is allowed 
to fall to very low values as suggested by simple flux tube rise 
simulations (case A), then the effect of diffusion is consider¬ 
able and the poloidal field is able to penetrate into the interior 
of the flux tube. On the other hand, if various effects in the top 
layers of SCZ keep the magnetic field above Beq (case Bl) and 
3Beq (case B3), then diffusion is less effective; in case B3 the 
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poloidal field remains confined in the sheath. For a poloidal 
field strength independent of depth, as assumed in these calcu¬ 
lations, the best agreement with the observed current helicity 
values is found for case B3. 

One rather strong prediction of our model is the existence 
of a ring of reverse current helicity on the periphery of active 
regions. On the other hand, the amplitude of the resulting twist 
(as measured by the mean current helicity in the inner parts 
of the active region) depends sensitively on the assumed struc¬ 
ture (diffuse vs. concentrated/intermittent) of the active region 
magnetic field right before its emergence, and on the assumed 
vertical profile of the poloidal field. Nevertheless, a mean twist 
comparable to the observations can result rather naturally in the 
model with perhaps the most plausible choice of assumptions 
(case B3). 

Thus, by studying the distribution of the azimuthal mag¬ 
netic field in sunspots and active regions, and by simulta¬ 
neously studying the subsurface magnetic field structure in 
emerging active regions by means of local helioseismology, it 
may be possible to test the present model and to throw some 
light on the conditions prevailing during the last phase of flux 
tube rise through the top layers of SCZ. We hope that mag¬ 
netogram data will improve significantly in the next few years 
and it will be possible to draw meaningful inferences. Since 
at present we have very poor understanding of the nature of 
the magnetic held or the effect of turbulence when a flux tube 
rises through layers immediately below the solar surface, indi¬ 
rect inferences from such observations of sunspots are of great 
importance. 

Lastly, we have used a simplifying assumption that the 
poloidal held that gets wrapped around the flux tube has a con¬ 
stant value of 1 G inside the SCZ. Dynamo models like the 
model of Choudhuri et al. (2005) can be used to calculate the 
distribution of the poloidal held inside the SCZ during different 
phases of the solar cycle. More detailed calculations of helic¬ 
ity can be done by using such poloidal held distributions in¬ 
stead of using the simple boundary condition (24). Since some 
groups have started reporting on the possible c ycle variation of 
helicity on the basis of th e observational data (IBao et all200(i 
iHagino & Sakuraill2005h . such calculations may become rele¬ 
vant in future when more detailed observational data are avail¬ 
able. 
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